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Abstract 

This paper derives and analyzes exact, nonlocal Langevin equations appropriate in a 
cosmological setting to describe the interaction of some collective degree of freedom with 
a surrounding "environment." Formally, these equations are much more general, involv- 
ing as they do a more or less arbitrary "system," characterized by some time-dependent 
potential, which is coupled via a nonlinear, time-dependent interaction to a "bath" of os- 
cillators with time-dependent frequencies. The analysis reveals that, even in a Markov 
limit, which can often be justified, the time dependences and nonlinearities can induce new 
and potentially significant effects, such as systematic and stochastic mass renormalizations 
and state-dependent "memory" functions, aside from the standard "friction" of a heuristic 
Langevin description. One specific example is discussed in detail, namely the case of an 
inflaton field, characterized by a Landau-Ginsburg potential, that is coupled quadratically 
to a bath of scalar "radiation." The principal conclusion derived from this example is that 
nonlinearities and time-dependent couplings do not preclude the possibility of deriving a 
fluctuation-dissipation theorem, and do not change the form of the late-time steady state 
solution for the system, but can significantly shorten the time scale for the approach towards 
the steady state. 
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I. Introduction 



Over the past two decades or so, a great deal of attention lias focused on the 
problem of understanding statistical field theory in a cosmological context, allowing 
correctly for the fact that the spacetime is not flat, or even static. This is a problem 
of interest both from an abstract fleld theoretic viewpoint and from a more practi- 
cal viewpoint which seeks to understand what the Universe was like at very early 
times. When addressing the latter sorts of issues, one is led oftentimes to implement 
various models and/or approximations which, although not justified rigorously, seem 
physically well motivated and even essential, if one wishes to obtain concrete results. 

One important ingredient in this sort of modeling is the idea that, at some level, 
the physical degrees of freedom of the Universe divide into two coupled pieces, a 
"system" component, the detailed evolution of which is for some reason of particular 
relevance, and a "bath" component, the detailed evolution of which is somehow ir- 
relevant. This sort of picture has arisen in at least four different settings. One such 
setting entails an understanding of inflation in terms of an inflaton fleld evolving in the 
rest of the Universe, which serves as an external environment or bath Another 
involves the general notion of "coarse-graining" as a physical mechanism in terms 
of which to extract quantum decoherence, this facilitating a "quantum-to-classical" 
transition in the early Universe Q. A third entails a more systematic development 
of statistical quantum fleld theory which uses a closed-time-path formalism to 
derive quantum dissipation and memory loss. And flnally, there is the intriguing, but 
not yet completely understood, program of stochastic inflation originally proposed by 
Starobinsky 0. 

Much work along these lines has been predicated upon the formulation of essen- 
tially ad hoc Fokker-Planck equations, which effectively introduce a bilinear interac- 
tion between the system and the bath. It seems crucial to understand the extent 
to which these sorts of heuristic models are in fact reasonable, i.e., approximately 
true in some appropriate limit, and, especially, how relaxing the basic assumption 
of bilinearity changes the underlying physics. In particular, is the standard sort of 
modeling legitimate if one allows, as in certain cases one must, for time-dependent 
couplings and frequencies and incorporates realistic nonlinearities? 

These are extremely difficult questions to answer in complete generality. However, 
as will be seen in this paper, they can be examined in certain cases by considering 
special models of systems coupled to baths comprised of time-dependent harmonic 
oscillators, where it is possible to derive exact, nonlocal Langevin equations simple 
enough to understand, both physically and mathematically. These models are nonlin- 
ear, time-dependent generalizations of phenomenological models which have proven 
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quite successful in other branches of physics, such as condensed matter physics or 
quantum optics. And, as in those settings, the models are well motivated phenomeno- 
logically, even if they are not derived ab initio. 

One of the objectives here is to derive exact generalized Langevin equations for 
these models, and in an appropriate Markov limit, their Fokker-Planck realizations. 
One is then poised to understand the sorts of new effects arising in these exact non- 
local equations which are absent both from more heuristic Langevin descriptions and 
from the exact nonlocal equations that can be derived for the special case of time- 
independent, bilinear couplings. What this entails is an analysis of the interplay 
between three different sorts of effects, namely linear and nonlinear noise, which may 
well have very different natural timescales, and the explicit time-dependence of the 
environment, reflecting the overall expansion of the Universe, which introduces yet 
another timescale. 

This analysis shows that (a) allowing for a nontrivial time-dependence necessarily 
induces qualitatively new effects like a mass (or frequency) "renormalization," even for 
the special case of bilinear couphngs 0); and, moreover, (b) allowing for nonlinearities 
in the system-environment coupling induces new effects aside from the usual "friction" 
term. One discovers, e.g., that nonlinearities give rise to an additional renormalization 
of the system potential, and that they imply a "memory" kernel which involves the 
state of the system. 

These results might suggest that, in the presence of such nonlinearities, one's naive 
intuition is completely lost. This, however, is not so: Even allowing for nonlinearities 
and time- dependent couplings in the interaction between the system and environment, 
one can, at least for the case of time-independent oscillators, where the bath may still 
be viewed as being "at equihbrium," still derive a simple fluctuation-dissipation the- 
orem 0. When the oscillators become time-dependent, the fluctuation-dissipation 
theorem will no longer be exact. However, this theorem does remain at least approx- 
imately true to the extent that the coupling between the system and environment 
is dominated by modes of sufficiently short wavelength. As will be shown below, in 
a cosmological setting this implies that, on scales short compared with the horizon 
length, it is still possible to speak of an approximate equilibrium and an approach 
towards that equilibrium. 

This paper focuses on obtaining a qualitative understanding of the effects of non- 
linear couplings and a time-dependent environment. A subsequent paper will present 
a concrete calculation, applying the technology of the time-dependent renormalization 
group to a simple cosmological phase transition. 

Section II of this paper focuses on the general problem of couplings between a 
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system and a bath, motivating in particular a rather general time-dependent Hamil- 
tonian which is amenable to a systematic analysis. Section III derives an exact, 
nonlocal Langevin equation from this Hamiltonian and then discusses its physical 
implications. All of this is completely general, not restricted in any way to a cosmo- 
logical context. Section IV then turns to a consideration of one specific cosmological 
model, deriving Langevin/Fokker- Planck equations for some collective degree of free- 
dom, such as the dilaton mode, evolving in a Landau- Ginsburg potential and coupled 
to scalar "radiation." This equation can provide one with a simple tool in terms of 
which to model a cosmological phase transition associated either with inflation or 
the formation of a cosmic string. Section V provides approximate solutions to this 
equation, which enable one to study the approach towards a (timc-dcpcndcnt) steady 
state. It is observed that such an approach towards "equilibrium" can be strongly 
influenced by the nonlinearities in the couplings ("multiplicative noise") and/or the 
effects of the time-dependent expansion of the Universe. 
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II. The System-Environment Splitting 

Two sorts of equations are ubiquitous in nonequilibrium statistical mechanics, 
namely (collisional) Boltzmann equations and Fokker-Planck equations. Boltzmann 
equations are appropriate for transport-type problems, involving strongly interacting 
particles, whereas Langevin/Fokker-Planck equations are useful in the study of phase 
transitions, Brownian motion, and the like. 

These two ubiquitous equations also find a place in cosmology. In the past, kinetic 
theory a la Boltzmann has been studied extensively, most prominently perhaps in the 
analysis of nucleosynthesis More recently, however, phase transitions have come 
into vogue, especially with regard to the infiationary scenario 0, and it is here that 
one encounters Fokker-Planck equations. 

That the standard methodology (with essentially trivial modifications) can be 
applied to the very early Universe involves a certain leap of faith; and while on the 
whole cosmologists appear comfortable with the status quo, doubts have certainly 



been voiced in the literature |jTO|]. The main objections relate to (1) the assumption 
of thermal equilibrium, (2) the absence of a clear separation of time scales, (3) the 
neglect of fiuctuations, (4) the validity of heuristic master equations, and, related to 
this, (5) the lack of a fundamental Liouville description, derived from a Hamiltonian. 
While it is fair to say that some of these objections have not been stated in a concrete, 
quantitative way, they do lead to a feeling of unease. 

In order to address some of these issues concretely, this paper will consider a proto- 
typical Hamiltonian which incorporates more or less realistic nonlinearities and time 
dependences, and then extract from that Hamiltonian an exact Langevin equation 
for the system variable. This is a Liouville approach, |10|, in which the fundamental 
equation is derived systematically without any ad hoc assumptions. However, the 
more difficult problem of justifying a full-blown nonlinear Boltzmann equation will 
not be treated here. 

The exact Langevin equation involves at least three distinct time scales: (1) the 
Hubble time tn, (2) the relaxation time t/j on which the system is affected by the 
surrounding environment, and (3) the time scale tc set by the decay of the noise 
autocorrelation function. If the system evolves under the infiuence of some nontrivial 
system potential U, there is also a fourth time scale ts, the time scale on which 
the system changes in response to U. The implementation of any approximation 
entails an assumption regarding the separation of these time scales. In particular, the 
possibility of an approximate local description, i.e., the existence of a Markov limit, 
depends critically on the assumption that tc ^ t/j, tn, and ts- An obvious point then 
is that nonlinearities, especially with respect to the system-bath coupling, can play 
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an important role by inducing multiplicative noise which can significantly alter t/j. 

The specific objective here is to motivate a phenomenological Hamiltonian for the 
system plus environment, and to analyze it rigorously to extract the physical effects 
contained therein. The system is taken to be some "collective coordinate" evolving 
under the influence of a "heat bath." Particular interest focuses on phase transitions, 
but the development here is in no way restricted to such a setting. Indeed, it should 
be stressed that this general approach has been used successfully in many other areas 
of physics, such as condensed matter physics, nonlinear optics, nuclear physics, etc. 

Given an arbitrary composite Hamiltonian, one wishes to introduce a splitting 
into a "system" piece, a "bath" piece, and an interaction term. If this split is to be 
useful, it must be true that in some sense the system is "small" compared with the 
environment. What this means is that, as far as the system is concerned, the full 
Hamiltonian 



where 5Hb is the Hamiltonian for a collection of harmonic oscillators and 5Hi an 



may well be one in which the oscillators are "fundamental" (e.g., the modes of some 
free fleld, as in black body photons), but this is by no means necessary: Assume that, 
in the absence of any coupling with the system, the environment is characterized by 
some flxed, possibly time-dependent, solution. Now allow for a weak coupling with 
the system, weak in the sense that each bath mode is only changed marginally. Then 
identify the g^'s as perturbed variables, i.e., degrees of freedom deflned relative to the 
flxed solution (e.g., phonons). This has two implications: (1) the environment can be 
visualized as a collection of oscillators with (possibly time-dependent) frequencies, so 
that 5Hb is quadratic in bath variables g^, and, (2) because the interaction of the 
environment with the system is assumed to be weak, in the sense that the individual 
bath modes are not altered signiflcantly, SHj must be linear in the qas. Note that 
one does not want to assume that the system is only weakly altered, so the interaction 
5Hi is not necessarily linear in the system variable x. In principle one can proceed 
without imposing any restrictions on the form of the system Hamiltonian Hs- 
Given the above set of assumptions, one can write that 



Ht = Hs + Hb + Hi 




is well approximated by 



H = Hs + 6Hb + 6Hi 



(2) 



interaction Hamiltonian linear in the oscillator variables qa [|n| |]T^. The heat bath 




(3) 




A 



(4) 
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and, in terms of relatively arbitrary functions F^, 



SHi = -J2Q^{t)TA{x,t)qA. (5) 

A 

It is, however, convenient to rewrite H in the manifestly positive form 

H=lv' + U{x, t) + W {pI + ^lit) [qa - Ta{x, t)f} , (6) 
where in terms of the "renormalized" potential Vren, 



f/(x,t) = Vren - - il^ it)T ^ (x, t) . (7) 
^ A 

Couplings of a system to some environment can induce finite and stochastic renor- 



malizations in the system potential, although this is not always so |T2[ . In this paper, 
the words "system potential" will always refer explicitly to the renormalized system 
potential. A physical restriction on the form of the couplings arises from the 
requirement that the renormalized potential not change the qualitative form of the 
bare potential. Thus, e.g., if the bare potential is a polynomial of order n, the renor- 
malization should induce no terms of order higher than n. This condition also insures 
stability of the system towards the destabilizing effects of multiplicative noise, since 
it implies that the stochastic forcing terms in the potential must be a polynomial of 
order < n. 



Finally, as emphasized, e.g., by Caldeira and Leggett |T^, it should be stressed 
that this is more than simply a toy model. This form of the Hamiltonian generally 
provides a correct description for any system which is only weakly coupled to its 
surroundings. To facilitate a concrete calculation, this Hamiltonian need only be 
supplemented by two inputs, namely the spectral distribution of the environmental 
modes and the form of the coupling to the system. For a general physical problem, 
these may either be extracted from experimental data or derived from theoretical 
considerations. 
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III. Time-Dependent Langevin Equations with Multiplicative Noise 

The equations of motion generated from the Hamiltonian @ clearly take the forms 



X = V, 

d d 

V = -^u{x,t) + J2Q^{t)[qA-rA{x,t)]—rA{x,t), 

U = Pa, 

PA = -n^[qA-TA{x,t)], (8) 

where an overdot denotes a time derivative d/dt. 

The fact that the equation for pA is linear in qa implies that one can immediately 
write down a formal solution for qA{t) in terms of r^(a;,t) at retarded times s < 
t. Indeed, let S'yi(t) and CAif) denote two linearly independent solutions to the 
homogeneous oscillator equation 

+ = 0, (9) 

chosen without loss of generality to satisfy the initial conditions C^(0) = S'a(O) = 1 
and S'yi(O) = Ca(0) = at some time t = 0. One then concludes exactly that 

QAit) = qA{Q)CA{t)+PA{Q)SA{t) 

+ /* ds n2{s)VA{x{s), S) [SA{t)CA{s) - SA{s)CA(t)] . (10) 

The integrand in ([T^) vanishes in the coincidence limit s ^ t. This, however, may 
be remedied by replacing Sa and Ca by —SaI^a ^"^^ —Ca/^a then integrating 
by parts. The net result is a formal solution 

qA{t)-TA{x{t),t) = [gA(0)-r^(x,0)]C^(t)+pA(0)^A(t) 

+ f dsWA{s,t)^TA{x{s),s), (11) 
JO as 

where the Wronskian 

WA{s,t) = SA{s)CA{t) - CA{s)SA{t). (12) 

Note that, for the special case of time-independent frequencies, C^(t) = cosVtAt and 

SA{t) = ^Ia^ COsilAt. 



By inserting (pT[) into the equation for i) and grouping terms suggestively, one 
then recovers an exact, nonlocal equation of the form 

'-T- ' 

A 



—Q^- I dsJ2Q2{t)WA{s,t)—rA{x{t),t)—rA{x{s),s) + Fs{t), (13) 



where 

d 



A 



{[g^(0) - Ta{x{0), 0)] CA{t) + PA{0)SAm . (14) 



In the spirit of the discussion in Section II, suppose now that the interaction 
between system and bath entails a polynomial coupling 

N 1 

r^(x,t) = E-7?^(^K, (15) 

where the functions 7^"* are arbitrary real functions of time. Equation (|13D then takes 
the form 

V = -— - / ds [K{t, s)v{s) + M{t, s)xis)] + Fs{t), (16) 
ox Jo 

where, in terms of the quantities 

AA{s,t) = (^^^^r\s)x-\s)^ (T.l^l\t)x-\t)^ (17) 

and 

BAis,t) = --^AAit,s), (18) 
n OS 

the "memory kernels" K{t, s) and M{t, s) are 

K{t,s) = Y.^^{t)AAis,t)WA{s,t), (19) 

A 

M{t,s) = Y.^^{t)BA{s,t)WA{s,t). (20) 

A 

The force Fs now reduces to 

F,(t) = j:^i{t)j:7^r\t)^"'-' {mo) - ta{x, o)] CA{t)+PA{o)SAm . (21) 

A rn 



Equation ([Tq) is considerably more complicated than an ordinary Langevin equa- 
tion. However, these additional complications need not preclude entirely the pos- 
sibility of a simple physical interpretation or the proof of a fluctuation-dissipation 
theorem. Provided that the oscillator frequencies Qa are not time-dependent, one 
can still prove a fluctuation-dissipation theorem, even if Ta is a nonlinear function 
of X 1^ and/or explicitly time-dependent. Indeed, consider an ensemble of initial 
conditions for which the first moments vanish identically, i.e., 

{QAm = (pa(0)) = 0, (22) 

with Qa{0) = Q'yi(O)— ryi(x(0), 0), and where the second moments are initially thermal, 
so that 

(pa(O)pb(O)) = nA{0)nB{0) (gA(O)QB(O)) = kBT5AB, (23) 

9 



where the angular brackets denote an initial ensemble average. One then computes 
exactly that 

{Fs{t)) = (24) 

and 

{Fsiti)Fsit2)) = kBTK{h,t2), (25) 

thereby identifying Fs{t) as a noise and proving a generalized fluctuation-dissipation 
theorem linking the noise autocorrelator with the "viscosity kernel" K{t, s). (Strictly 
speaking K{t, s) need give rise to a true viscosity only in the case of an ohmic envi- 
ronment, a point which will be discussed later.) Equation ( [TE| ) can now be viewed as 
a nonlinear, nonlocal Langevin equation. 

That a fluctuation-dissipation theorem can hold even in these more complicated 
settings is a formal consequence of the fact that the nonlinearities and time- dependences 
in (p!6|) enter into the memory kernel K{ti, and the autocorrelator {Fs{ti)Fs{t2)) in 
exactly the same way. Physically, this result can be understood as follows: If the basic 
picture is valid, the total energy is dominated by the constant energy of the heat bath 
(recall that one is now assuming that the bath frequencies are time- independent), so 
that energy is approximately conserved, even if the system Hamiltonian Hs is time- 
dependent. Since fluctuations induce a monotonic increase in the system energy, there 
must be some source of dissipation if that energy is to remain bounded. However, 
if that dissipation is too strong and dominates the fluctuations, the system energy 
will vanish at late times, which is clearly unphysical for a system coupled to a finite 
temperature heat bath. The fluctuations and the dissipation must clearly balance if 
the system is to have a finite but nonzero energy at late times. 

It should be observed that the Langevin equation (^) reduces to a well known 
form in an appropriate limit. If one neglects all nonlinearities in the coupling between 
system and bath, assuming that Va oc x, one immediately recovers a special model 
considered previously And, moreover, if one assumes further that VLa and Va 
are independent of time, one is reduced to the well known independent oscillator 
model [|^. It is thus possible to address systematically the question of how the 
incorporation of nonlinearities and/or time-dependences leads to systematic changes 
in the Langevin equation derived for that original model. 

When one neglects both the nonlinearities and the time dependences, the mem- 
ory kernel M(s, t) vanishes identically and, moreover, the remaining memory kernel 
K{s, t) contains no explicit x-dependence. The stochastic force Fg involves x only 
linearly, through the propagation of an initial condition. It thus follows that one 
recovers a relatively simple equation involving a (nonlocal) friction oc v{s) and purely 
additive noise. To the extent that the function K{t, s) is sufficiently sharply peaked 
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about the coincidence limit s — > t (what precisely this entails will be discussed below) , 
one can then approximate v{s) by its value at time t, in which case one obtains a 
Markovian equation involving a friction oc v{t). 

When one allows for a nontrivial time- dependence in the oscillator frequencies or 
the coupling between the system and bath, but as yet no nonlinearities, an additional 
nonvanishing kernel M{s, t) appears. In a Markov approximation, this leads to a 
new term in the Langevin equation proportional to x{t) which corresponds to a time- 
dependent change in the system potential. In a field theoretic context, this would 
be interpreted as a mass renormalization. In this case, Fg still gives rise to additive 
noise, but that noise acquires an explicit time-dependence. 

When instead one incorporates nonlinearities but no time-dependences, the mem- 
ory kernel M{s, t) still vanishes, but the other kernel K{s, t) becomes significantly 
more complicated, involving not simply an autocorrelator for the mode functions 
Ca oc cosQAt, but a correlator of dTA/dx with itself. In other words, the nonlinear- 
ity implies that the evolution of x actually involves an x-dependent memory. This is 
hardly surprising. Indeed, as one might have anticipated, e.g., by analogy with the 
theory of the Brownian motion, the integral is nothing other than the autocorrela- 
tion function for the forces associated with the interaction of the system with each of 
the bath modes. In this case, Fg also acquires an additional a;-dependence implying 
that the noise will depend not only on the initial conditions x(0), ^^(0), and Pa(0), 
but upon the x-dependent state of the system as well. In other words, the noise is 
multiplicative. 

It should be observed that the memory kernel M{t,s) does not enter into the 
fluctuation-dissipation theorem for a time-independent bath. Its only effect is to 
induce a systematic renormalization of the system potential, which will of course 
affect the form of the late time solution. 

One further point should be stressed. The fluctuation-dissipation theorem of equa- 
tion (|25|) refers explicitly only to the total force Fg and the total memory kernel 
K{t,s). However, it is easy to see that analogous theorems also hold separately for 
each term oc 7^''(s)7^^(t). In this sense, the fluctuation-dissipation theorem is truly 
microscopic. 

Equation (^) is an exact, nonlocal equation. Only to the extent that this equation 
can be approximated as Markovian can one derive from it a Fokker-Planck equation. 
Such a Markov limit implies (1) that the memory kernels K{t, s) and M{t, s) may be 
approximated as essentially local in time, and (2) that the autocorrelation function 
for Fg falls off rapidly as |t — s\ increases. It is clear by inspection that, if K{t, s) is 
essentially local, so is M{t, s). Further, the fluctuation-dissipation theorem guarantees 
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that, if K{t, s) is local, the noise autocorrelation will be as well. In considering the 
validity of a Markov approximation, it thus suffices to consider K{t, s). 

The Markov limit can oftentimes be justified approximately for smooth spectral 
distributions, given an appropriate separation of time scales. Let tc denote a char- 
acteristic time scale on which the memory kernels K{t, s) and M(t, s) decay, and let 
tgys denote a characteristic time scale on which the system velocity v and/or position 
X change significantly, either in response to the environment (t/j) or to the system 
potential U{ts). (Provided that the mode functions are all oscillatory, K{t, s) will de- 
cay much faster than |t — s| at late times, so that one can clearly identify a time scale 
tc-) In a time-independent setting, the Markov limit then follows when tc -C tsys- 
If the oscillators are time-dependent, there is another relevant time scale, tn, the 
characteristic time on which the frequencies change. To the extent that tn is much 
larger than both tgys and tc, one anticipates that the time- dependence may be viewed 
as a perturbation, and that the Markov limit should still obtain. If tc -C tgys, but 
is not short compared with t^, the time dependence can no longer be viewed as a 
perturbation, but it may still be true that a Markov approximation can be justified. 
This should, e.g., be the case if all the modes still oscillate and/or the coupling of 
the system to the longest wavelength modes is relatively weak. If, however, tn is 
not much larger than tc, one expects that a local Fokker-Planck description will be 
inappropriate. This is, for example, true for the specific example discussed in Sections 
IV and V. 

For the special case of "ohmic" environments, the nonlocal term in the exact 
Langevin equation involving K{t, s) reduces to the usual linear viscosity seen in 
heuristic ad hoc Langevin descriptions. These environments are characterized by 
time-independent frequencies and have a spectral distribution g{fl) oc Q"^, with an 
upper cut off at some Qmax, and all the oscillators are assumed to couple bilinearly 
to the system with an equal strength [|12| ||^. If the coupling is not bilinear, one still 
recovers a viscosity that is linear in v, but this viscosity will dependent explicitly on 
X. If the spectral distribution g{Q) differs only slightly from oc Q"^, that difference 
may be treated perturbatively to extract calculable modifications in the form of the 
Markovian equation. However, for spectral distributions which are very different, e.g., 
"supra-ohmic" distributions oc Q'^, the local Langevin equation can be higher order 
in time derivatives and need not contain a simple viscosity term oc f (t). One concrete 
example thereof is provided by an electron interacting with its self-electromagnetic 



field [0] . However, for the cosmological example considered in Sections IV and V, the 
spectral distribution will be nearly "ohmic," so that the local Langevin description 
will contain an ordinary viscosity and thus admit a Fokker-Planck realization. 
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IV. A Cosmological Example 

The objective of this Section is to formulate a nonlocal Langevin equation in terms 
of which to describe the evolution of some system variable, evolving in a Landau- 
Ginsburg potential V^en and coupled to a bath of scalar black body "radiation." The 
aim of this Langevin equation is to provide a quasi-realistic model for a cosmological 
phase transition. The entire analysis will be classical, allowing for thermal fluctuations 
but not for quantum fluctuations. This may prove appropriate either in the context 
of some versions of inflation or, alternatively, in the formation of cosmic strings or 
baryogenesis. Attention here focuses on formulating the problem and discussing its 
physical potentialities. Quantitative details will be provided in a subsequent paper. 

In what follows the bath will be idealized as a collection of oscillators, character- 
ized by time-dependent frequencies appropriate for a scalar field with a general C,R 
curvature coupling. The case = 1/6 corresponds to conformal coupling, whereas 
^ = yields minimal coupling. Consistent with the discussion in Section II, the 
interaction Hamiltonian will be taken as linear in the oscillator variables qa, but it 
can involve an arbitrary quadratic dependence on the system variable x, with both 
linear and nonlinear pieces. By allowing for both linear and nonlinear couplings, and 
varying the relative strengths of these two different contributions, one will be able 
to compare the effects of additive and multiplicative noise. The analysis will be ef- 
fected in the conformal frame, in terms of a conformal time coordinate rj satisfying 
dr] = a~^dt, where a denotes the scale factor. It will, moreover, be assumed that 
the spatial curvature of the t = constant slices vanishes, so that one is considering a 
k = Friedmann cosmology. 

Given these assumptions, one is led directly to a Hamiltonian of the form moti- 
vated in Section II, namely 

H = ^v' + U{x, V) + 1^ {pI + ^liv) [Qa - TaIx, v)]}' , (26) 

where, in terms of constants A, 6, and a, 

f/(x, v) = \mx' + l^ivW' + ^^(^7)x^ (27) 

and the coupling 

r^(x,r7) =7«(r7)x + 7?^(r7)x2. (28) 

Recall from equation (0) that U and the Landau- Ginsburg potential Ken are con- 
nected by a term involving the couplings Ta- It follows that, if the 7^^^'s and 7a^^'''s 
are both nonvanishing, this coupling will in general induce a cubic term in U, al- 
though the coefficient 6 of that term could vanish. The coefficients A and a in U are 
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also "dressed" quantities, differing from the bare quantities in V^gn because of the 
terms involving the F^'s. At this stage, the explicit time dependence of A and a and 
the couplings 7^^^ and 7^"* may be treated as more or less arbitrary. The frequencies 
satisfy 



a" 



-(1-60- (29) 

where a prime ' denotes a conformal time derivative dldr]. 

Suppose now that, in terms of cosmic time the scale factor a evidences a simple 
power law time dependence a — t^, with p > 1/2. It then follows that 

a(x{r)-rjoy/^^-P\ (30) 

where 770 denotes an integration constant, so that 

where 

is intrinsically positive when p > 1/2 and 1 — 6^ > 0. The quantities lo^ denote 
eigenvalues of the spatial Laplacian. 

To proceed further, one needs to determine the mode functions for the time- 
dependent frequencies. These can clearly be evaluated in terms of Bessel functions. 
However, in so doing there are at least two possible ways in which to proceed, namely 
considering (1) complex modes involving Hankel functions or (2) real modes involving 
ordinary Bessel and Neumann functions. This paper will adopt the (less conventional) 
second choice, since it provides for a more direct connection with earlier work on 
Langevin equations: in the absence of a time-dependent expansion, the mode func- 
tions reduce to sines and cosines. Straightforward algebra reveals that the equation 

,2 

Ea^O (33) 



■—4 ^ 



2 



is solved by a general 

5a = [i^A iv - Vo)]^^'^ {uja {v -Vo)), (34) 
where denotes an arbitrary solution to Bessel's equation of order 

"^-(l-^of^-I-'-^ + i. (35) 

Note that, for the special cases = 1/6 (conformal coupling) and/or p = 1/2 
(a Universe dominated by conformal electromagnetic radiation), = 1/2 and the 
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solutions reduce to sines and cosines. For these particular values, the bath Hamil- 
tonian 5Hb is time-independent in the conformal frame. This implies that the physical 
frequencies are simply red-shifted uniformly as the Universe expands. 

Presuming that the initial value problem for the Langevin equation is specified at 
time ?7 = 0, the appropriate solutions will be Ca and Sa, combinations of Bessel 
and Neumann functions Jy and A^,^, modulated by factors ujA{f] — ^Q)i satisfying 
^^(0) = dSA{Q)/dr] = 1 and Sa{0) = dCA{0)/dr] = 0. If the functions 

Ca and Sa are, for long wavelengths with \u!a ij] — Voll ^ Z^; 'vsry different from sines 
and cosines. Indeed, they are not even oscillatory. However, to the extent that such 
longer wavelengths do not couple significantly to the system, one can approximate 
the Sa's by the forms appropriate when \uja ij] — ?7o)| ^ In this case, the mode 
functions reduce to 



'a(^) ~ sin UA iv -Vo)-^(^f^+ ^) 



• 2x1/2 

-a(?7) ~ ( - cos 



TT 



2 1^+2 

TT / 1^ 



(36) 
(37) 



L^A iv -Vo)-^[^^+ 2) 
i.e., ordinary sines and cosines, modulated by phase shifts which can of course be 
absorbed in the normalizations. 

It is clear that, in this limit, one recovers a relatively simple nonlocal Langevin 
equation, for which the only explicit time dependences are in the potential U and 
the couplings between the system and the environment. It thus follows that, in this 
approximation, a fluctuation-dissipation theorem holds, so that one would antici- 
pate an evolution towards some steady state solution at late times. The fluctuation- 
dissipation theorem is of course exact when ,^ = 1/6 and/or p = 1/2. 

It should be stressed that the condition \uja{v — ^70)! ^ A* has a very simple physi- 
cal interpretation. Reexpressed in terms of the physical cosmic time t, this condition 
becomes {ujA/ci)t ^ — p\, this corresponding, for /i and p of order unity, to the 
demand that the physical period of the oscillation be short compared with the time 
scale tn associated with the expansion of the Universe. In other words, the wavelength 
must be short compared with the horizon length. 

In general, however, a fluctuation-dissipation theorem does not hold. Recall that 
the memory kernel K{t, s) satisfies 

K{t,s) = Y.Q^{t)AA{s,t)WA{s,t), (38) 



in terms of the Wronskian of equation ([T^). Alternatively, a thermal average of the 
noise autocorrelation, taken at the initial time rj = 0, satisfies 

(kBT)'' {F,{t)Fs{s)) = J2n^{t)AA{s,t)QA{s,t), (39) 

A 
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where 



+ SA{t)SA{s) 



(40) 



When the frequencies Qa are all time-independent = 1/6 or p = 1/2), VFA(s,t) 
and Qa{s, t) are in fact equal, so that one recovers the fluctuation-dissipation theorem 
(p5|). In general, however, this is clearly not so; and the typical size of the fractional 
deviation between these two quantities provides a concrete measure of the degree to 
which the fluctuation-dissipation theorem fails. Suppose that Ucr denotes a typical 
frequency associated with the coupling of the system and the environment. It then 
follows straightforwardly that the fractional amplitude of the deviation is of order 

7 7 7- (41) 

The fractional deviation from a fluctuation-dissipation theorem scales as the ratio of 
a characteristic oscillator period ~ {ucr/a)'^ to the expansion time scale tn- 

It remains to consider the circumstances under which a Markov approximation 
can be justifled. Note flrst of all that, to the extent that the bath oscillators are in- 
terpreted as representing the modes of some free fleld, one would anticipate a spectral 
distribution oc Ua, so that, given some cutoff ujmax, one can pass to a continuum limit 



Y,-^ dco iu^. (42) 

Consistent with the equation for 5Hi, now separate out the explicit ^-dependence in 
the system-bath couplings and write 7^"^ (77) = c^2\v) I^a{v)- 

As a particularly simple flrst approximation, suppose that c^"* is essentially in- 
dependent of frequency. And, moreover, assume that the longest wavelength modes 
are not very important in the coupling, so that one can neglect the frequency shift 
associated with the expansion of the Universe and set VLa ~ ^a- Given these approx- 
imations, one has, e.g., that 

K{r],s)^ f dujJ2c^''\s)c^""\v)x'"\s)x"''\v)cosuj{r]- s). (43) 

m,n 

Suppose, however, that one can also neglect the time-dependence of the c's, a reason- 
able assumption, e.g., if the c's change only on an expansion time scale iff. One can 
then effect the du integration explicitly to obtain (for large uJmax) a delta function 
^d{v It follows that the memory of K{t, s) is indeed very short, so that, presum- 
ing that the time scale ts associated with U is not too short, one can approximate 

1 2 



K{v,s) 



5o{r,~s)=21C{x,r,)5D{r,-s). (44) 



16 



An analogous expression holds for the other memory kernel M{7], s). 

If the ca's depend strongly on frequency and/or the spectral distribution is con- 
siderably different, the analysis becomes more complicated. It is, however, evident 
that, whenever the couplings are such that the lowest frequency modes are not very 
important, so that tc ^ and tn, and the distribution of modes is approximately 
ohmic, a Markov approximation should in fact be justified, leading to a viscosity 
oc JCv. 

Under these circumstances, the nonlocal Langevin equation can be well approxi- 
mated by a local equation of the form 

dU 



V = - M{r])x{ri) - }C{r])v{r]) + Fs{r]), 



where, for an initial thermal ensemble, 

{F,{7])F,{s)) = 2kBTIC{x,7])6D{^ - s) 



(45) 



(46) 



This local Langevin equation leads immediately to a Fokker-Planck equation of the 
form „ , , ^ 

-— - Mx- JCv] f -fcoTr ^ 
ox 



dr] dx dv 



0. 



(47) 



Translated back into the physical frame, the Langevin equation (^5]) becomes 



V 



1 dU 



M{t)X 



3- + K(t) 

. a 



a^dX 

where the overdot denotes differentiation with respect to cosmic time, and 

X 

a 



X 
V 

K(t) = 

m{t) = 

F.(t) = 

The noise autocorrelator is now 



if - — - — 



a 

M 

a? 

El 



a" 



a" 



- +-, 



(48) 

(49) 
(50) 
(51) 
(52) 
(53) 

(54) 



Note that, in the physical frame, there are two sources of damping, namely the 
viscosity oc K.V and the cosmological frame- dragging cc H = a/a. The corresponding 
Fokker-Planck equation is now 



d 



1 dU 
'^dX 



M{t)X 



K(t) V 



0. 



(55) 
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The form of this equation is consistent with the Fokker-Planck equation analyzed 
by Brandenberger et al [|l|. Given, however, that the approach followed here is entirely 
different from that of Ref. [|I| a more detailed comparison is appropriate: The present 
formalism is intended primarily to describe the behavior of a homogeneous degree 
of freedom (dilaton mode). However, for the special case of a free field with no 
mode couplings, the Fokker-Planck equation (^) holds equally well for any single 
field mode coupled to a heat bath, which is precisely the model problem considered 
by Brandenberger et al. Comparing (|55|) with equation (20) of Ref. one finds 
that, even though the latter equation was arrived upon in a rather different way, the 
diffusion terms are consistent in that they scale the same way with a, and that in 
both cases the "Hubble damping" term is present. However there is one important 
difference. While the present model incorporates back-reaction effects due to the heat 
bath, namely the normal viscosity K and the potential renormalization M, these 
physical effects are not taken into account in Ref. A quantitative assessment of 
the importance of these effects will be presented elsewhere. 

Having formulated the Fokker-Planck equation, it is worth recalling once again 
the critical assumptions that went into its derivation. (1) The decay time tc must 
be short compared with the time scales and ts on which the system changes in 
response either to its potential U or in response to the environment and with the 
expansion time tn- (2) For 00(1] — tjq) <^ /i, the modes are non-oscillatory. One must 
also assume that the coupling of the system with these infrared modes is negligible, 
so that there is no significant long time tail to K{ri, s) to prevent the existence of 
a Markov limit. (3) The spectral distribution must be approximately "ohmic," with 
g{uj) oc cu^, so that the nonlocal contribution involving Kit], s)v{s) gives rise to an 
ordinary viscosity. Fortunately, this is precisely what one expects of a heat bath of 
thermal photons. 
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V. An Approximate Solution 



Given the assumed form (|26D for the Hamiltonian (with the system-bath coupling 



assumed to be weakly time-dependent), the coefficients /C and M. may be written as 

/C = Ao + 2AiX + Aax^ (56) 



and 



Ai = fiQ + 2fiix + ii2X^ (57) 



and the noise correlator 

(F,(r/)F,(s)) = 2kBT}CSDiv ' s) = 2kBT [Aq + 2AiX + Aax^] 5d{v - s), (58) 

where the A„'s and independent of x and v. Note that Aq and A2 are neces- 

sarily positive definite, but that Ai and the of indeterminate sign. /C gives 

rise to an x-dependent viscosity, whereas Ai induces a further renormalization of the 
potential U in terms of new coefficients A, B, and S: 

1 1 1 12 1 

^ i^, V) = 2^{v)x'^ + ■^^{v)^^ + 2 ^(^)^^ + 2^ox'^ ^ ^^^^^ ^ i'"^^^ 

= ^A{r])x^ + ^e{7])x^ + ^J:{r])x'^ (59) 



By inserting (^) and ( pSf ) into (|47|) , one obtains an explicit Fokker-Planck equa- 
tion, which one may hope to solve. Unfortunately, however, in general such a multi- 
variable Fokker-Planck equation cannot be solved exactly except via numerical tech- 
niques. Nevertheless, one can at least obtain an approximate solution for the ex- 
pectation value (E) of the system energy E using the so-called "energy envelope" 



technique introduced by Stratonovich [|T^] and further developed by Lindenberg and 
Seshadri ||T6| . 

The idea underlying this approximation is in fact straightforward: In many cases 
of physical interest, such as that considered here, one can visualize the system as 
exhibiting (nonlinear) oscillations on the time scale ts ~ ujq'^ associated with the 
potential U, oscillations which are eventually altered by the coupling with the envi- 
ronment on the damping time scale tR. To the extent that the damping time and 
the expansion time tn are both long compared with ts, one can then assume that the 
system energy E is nearly conserved on time scales ~ ts, and varies only on a time 
scale much longer than the time associated with variations in x. It is thus natural to 
transform from x and v to new variables x and E, to treat E as an adiabatic invariant, 
and to implement an "orbit average" of the transformed E-x Fokker-Planck equation 
to extract an equation involving only E and r]. 
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Implement, therefore, a change of variables from (x, v) to (x, -E), where, explicitly. 
E = jv'^ +U{x,ri), to obtain a new Fokker-Planck equation for W{x,E,t) satisfying 



/(x, f , rj)dxdv = W{x, E, rj)dxdE. 



(60) 



To the extent that the energy is approximately conserved during a single oscillation 
of the system, one can then assume that 



W{x,E,r]) = {2^'{E) [E -l{{x,r])]^^^}-\Wi{E,r]) 



(61) 



where 

$(E) = Jdx[E- U{x, r/)]^/^ (62) 

and ' now denotes a d/dE derivative. Here the integration extends over the values 
of X along the unperturbed orbit associated with E. Note that the prefactor of Wi is 
simply the relative amount of time that, for fixed E, the system spends at each point 

X. 

By integrating the Fokker-Planck equation for W{x, E, rj) over x, one obtains the 
desired equation for Wi{E,ri), which takes the form 



-(- 

\dE 



X'{E) 

+X2{^{E)-kBT^'iE))}] 



{Ao mE) - kBT^\E)) + 2Ai {x{E) - kBT^{E)) 



+kBT 



dE^ 



where 



and 



X'{E) 
xiE) 



{Xo^E) + 2\ixiE) + \2^{E)} 



W,{E,r^), 



(63) 



ldxx[E-Uix,,)] 



1/2 



11/2 



(64) 
(65) 



^(E) = J dxx^ [E-U{x,r])f 

Note that this differs from equation (3.8) in ||16|, which assumes implicitly that the 
system's unperturbed orbit is symmetric about x = 0, so that (0) vanishes identically. 

Unfortunately, for a generic potential U the functions $(-E), x{E), and "^{E) 
cannot be evaluated analytically, so that one cannot realize the right hand side of 
(|53|) explicitly in terms of simple functions of E. Thus, e.g., for the quartic poten- 
tial (|60D, these functions can only be expressed as elliptic integrals, which must be 
evaluated numerically. There is, however, one limit in which one can proceed ana- 
lytically, namely when the energy E is small and the system is oscillating about a 
local minimum of U. The obvious point is that, in this limit, one can evaluate the 
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orbit integrals, assuming that the system is effectively evolving in a simple harmonic 
oscillator potential. 

Suppose for simplicity that G = 0, so that the dressed potential lA is itself of the 
Landau-Ginsburg form. When S is positive, one can then approximate the system as 
oscillating with squared frequency VLq = S about the origin. And similarly, when S is 
negative, and W is a "Mexican hat" potential, the system can be assumed to oscillate 
with squared frequency Qq^ = 2 |S| about one of the two minima at xq = ±(A/ 

In either case, one can evaluate $, x? ^''^^ ^) to realize the right hand side of 
(|53|) in terms of polynomials at most quadratic in E. And, given such an explicit 
representation, it is straightforward to derive from the Fokker-Planck equation a 
transport equation involving the time derivative of the first energy moment 

{E{r])) = JdE E Wi{E,7]). (66) 

Suppose in the first instance that S is positive, and that the system is executing 
small oscillations about the ground state x = with uJq^ = S. In this approximation, 
one verifies that, to the extent that the time dependence of E, ujq, and the A„'s may 
be ignored, the moment equation takes the form 

I- {E{r^)) = ksTXo - f Ao - {E{r^)) - ^ (E\r^)) . (67) 



OT] \ UJq J CUi 







Note that, because of the reflection symmetry x — > —x for the potential U, the 
functions x = x' = O5 so that the contributions involving Ai vanish identically. 

Unfortunately, this equation still cannot be solved exactly for {E{ri)), as it involves 
the unknown function {E'^{ri)). To obtain a formula for {E'^ij])), one must consider 
the second moment equation, which in turn relates d {E'^{ri)) /drj to the third moment 
{E^{ri)). In the spirit of (say) the BBGKY hierarchy, one requires a truncation 
approximation. 

As in Refs. |]15| and |jl6|, suppose therefore that 

V(t)) ^ K{E{t))\ (68) 



with K = 2. One knows that, when the system is at equilibrium, with energy E = ksT, 
this equation is satisfied identically for k = 2, and one might expect on physical 
grounds that, before the system is "at equilibrium," the energy distribution will be 



narrower and k < 2. As emphasized by Lindenberg and Seshadri [0, this truncation 
approximation thus leads to an upper limit on the time scale on which the system 
"equilibrates" with the bath. Given this truncation, one can immediately write down 
the solution [|1^] 

^ ksT (Eq + AkBT) - AkT {ksT - Ep) exp {- [{A + 1)/A] XqT]} 
^ ^"^'^ iEo + AkBT)~ikBT-Eo)exp{-[iA + l)/A]\or^} ' ^^^^ 
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where 

(™) 

In the hmit that A2 0, the muhiphcative noise "turns off" and the system 
approaches an "equihbrium" with {E) = ksT on a time scale tn ~ Aq ^. If A2 7^ 0, 
the system still evolves towards an equilibrium with [E) = ksT, but the time scale 
tu can be altered significantly. Indeed, in the limit that Aq 0, the additive noise 
"turns off" and tn ~ uJq^/ {\2kBT). It thus follows that, when the nonlinear coupling is 
sufficiently strong, the system may be driven towards equilibrium, not by the ordinary 
additive noise associated with the linear coupling, but primarily by the multiplicative 
noise associated with the nonlinear coupling. 

The only point that remains to be checked is that one is still assuming, as is 
implicit in this "envelope" approximation, that the time scale uJq^ is much shorter 
than the damping time. This, however, is clearly the case when Aq and A2 are not 
too large. Indeed, one verifies that (a) the weak damping approximation is legitimate 
but (b) multiplicative noise dominates the evolution towards an equilibrium whenever 

m 

Turn now to the case when S < and the system is oscillating about one of 
the two minima of the potential xqt^O. This is the case relevant to first order phase 
transitions. Here the terms involving x and its energy derivative do not vanish, and 
the formulae for $ and \I' acquire additional terms involving the location xq of the 
new minimum. However, one still recovers a relatively simple exact equation for 
d {E) /drj. Specifically, one finds that, in this case, ( p7| ) is replaced by 

I (E(,)> ^ tsTL - - ^) iSM) - ^ (S^,,)) . ,72) 

where now 

S = E-U{xo) (73) 
denotes the system energy defined relative to the minimum of the potential, and 

L = Ao + Aixo + AqXq^ (74) 

plays the role of a "dressed" Aq. The obvious point here is that, since one is effectively 
expanding in a Taylor series around the point xq = ±(A/ |S|)^/^, the coupling terms 
oc Ai and A2 will induce ( X — Xq )-independent effects. 

Note in particular that, if |xo| is large, as will be the case when S <^ A, the dressed 
L can be much larger than Aq. This implies that, in this case, the nonlinear couplings 
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can reduce the overall equilibration time both through the introduction of a new term 
oc A2 {E"^) and through an increase in the effective linear coupling. 

In any event, to the extent that U{xo, rj) is only slowly varying in time, the deriva- 
tive d {E{ri)) /drj can be replaced by d {S{r])) /drj. And, to the extent that the co- 
efficients L and A2 may be approximated as time-independent, (^) can be solved 
analytically. The result is an expression identical to (p3|), except that E is replaced 
by the shifted £ = E — Uq and 

If the energy E and the couplings cannot be treated as independent of time, 
the analysis becomes more complicated, but, at least when E is small, one can again 
formulate an analogue of (^). If E depends explicitly on r], the moment equation 
will of course acquire an additional term {BE /drj). Suppose that the system may be 
approximated as simply oscillating with squared frequency T,{rj) about Xq = 0. One 
then concludes that 

Consistent, however, with the Ansatz (|61|), can be replaced by its "orbit aver- 

aged" value E/2 (this is the expectation value associated with the W of eq. (|6TD), so 
that the moment equation may be written in the form 

d 1 . „ . „, kBT\2\ A2 /^2. 



^ {E{V)) = 2 (E{v)) ^InS + ksTXo - ^Xo - j {E{v)) - ^ {E^v)) ■ (77) 

Suppose now that the time dependence of S is relatively unimportant, i.e., that 
the time scale on which S changes is long compared with the time scale on which 
the environment effects the system. In this case, it makes sense to speak of (at 
least) an approximate approach towards equilibrium on a time scale set by the time- 
dependent couplings. If the coupling between the system and the environment is 
dominated by Aq, one thus infers a decay of initial conditions and an approach towards 
an equilibrium, driven by the additive noise, proceeding as exp[— / drj/ Xq^t])]. And 
similarly, if the coupling is dominated by A2, one has an approach towards equilibrium, 
driven by the multiplicative noise, proceeding as exp{— / dr][T,/ X2{kBTri)]}. 

Suppose, however, that the time dependence of S is important, and that S changes 
appreciably on time scales -C t^. In this case, one can no longer speak of a simple 
approach towards equilibrium, since the form of the system Hamiltonian is actually 
changing on a time scale -C t_R. One now concludes that, in a first approximation, 
E(ri)/E(0) = [E(77)/E(0)]^/^, and that the coupling with the environment is only a 
perturbation on this simple power law evolution. 
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VI. Conclusion 



The principal thrusts of this paper have been (a) the motivation of a relatively 
general phenomenological Hamiltonian, in terms of which to characterize the evolution 
of a "system" degree of freedom with its surrounding environment; (b) the rigorous 
derivation of Langevin and Fokker-Planck equations for a system described by this 
model Hamiltonian; and (c) a qualitative analysis of the new effects incorporated in 
these equations which are absent from other, more heuristic, descriptions. The model 
considered in this paper has several potential applications, the most obvious being 
to the study of cosmological phase transitions. Unfortunately, however, the Langevin 
and Fokker-Planck equations derived here are, except in a few simple cases, very 
difficult to solve analytically. For this reason, a numerical study applying the model 
to several cosmological problems of interest, such as the onset of new inflation and 
first order cosmological phase transitions, is currently underway. 
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